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CONVEX:  A  COMPUTER  PROGRAM  FOR  SOLVING  CONVEX  PROGRAMS 

by 

H.  0.  Hartley,  R.  R.  Hocking,  L.  R.  LaMotte,  H.  H.  Oxspring 

Introduction 

The  purpose  of  this  technical  report  is  to  describe  a  computer  program  imple¬ 
menting  the  convex  programming  algorithm  developed  by  Hartley  sad  Hooking  (1963). 
The  program  has  been  used  successfully  on  a  number  of  problems  connected  with 
Project  Themis  as  well  as  numerous  other  problems  which  h8ve  arisen  at  Texas  A&U. 
The  program  is  also  being  used  by  a  number  of  other  institutions  as  shown  on  the 
attached  List  of  Users. 

The  report  is  in  two  parts.  Part  I  is  a  description  of  the  algorithm  as 
extracted  from  the  initial  paper  and  Part  II  gives  a  documentation  of  the  program. 
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PART  I:  Description  of  the  Algorithm 
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INTRODUCTION 

The  general  mathematical  programming  problem  can  be  described  as 
maximize  g(x) 

subject  to  f^(x)  <  0  i»l  . ..  n 

The  convex  programming  problem  is  an  important  special  ease  which  is  amenable 
to  mathematical  and  numerical  treatment  and  requires  that  g(x)  is  concave 
and  the  f^(x)  are  convex,  real-valued  functions  of  the  n-vector  X  for  all 
real  X.  The  non-negative  condition  on  the  variables  may  be  included  but 
it  introduces  no  loss  of  generality.  It  is  also  assumed  that  the  functions  ere 
differentiable.  Certain  special  problems  exist  in  which  these  properties  hold 
only  in  a  subset  of  the  n-space  but  they  can  frequently  be  handled  by  slight 
modifications  of  the  algorithm  to  be  discussed  in  the  following. 

Formulation  of  the  Convex  Problem  as  a_  Linear  Problem 

The  problem  given  above  is  equivalent  to  that  of  maximizing  the  co-ordinate 
Xn+1  subject  to  the  originel  restrictions  and  the  additional  restriction 
Xn+1  -  g(x)  <  0.  Note  that  the  concavity  of  g(x)  implies  that  this  new 
restriction  is  convex  and  oould  just  be  adjoined  to  the  original  set  as 
fm+i(x)  <  0.  It  does  have  particular  significance  however  and  so  it  will  be 
carried  in  its  original  form. 

This  new  problem  has  a  convex  feasible  region  in  (n+l)  -  space  bounded 
on  the  sides  by  cylindrical  Bets  formed  from  the  boundaries  f^(x)=0  and  x=0 
and  bounded  from  above  by  the  surface  xR+^  =»  g(x) .  This  feasible  region  is 
now  approximated  by  a  new  one  which  is  defined  by  hyperplanes  in  (n+l)  -  space, 


~y 


that  ig ,  the  convex  feasible  region  defined  by  the  restrictions 
f^(x)  ^  0  i  «  1  ...  ffl  +  1  is  replaced  by  one  vhicii  is  defined  by  chs 
intersection  of  a  large  set  of  linear  half-spaces.  These  hyperplanes  are 
constructed  as  follows : 

(i)  Impose  an  arbitrarily  fine  grid  on  n-space. 

(ii)  Construct  a  tangent  plane  to  each  surface  x  +  ^  «  f^(x) 
i  »  1  ...  a  +  1  at  each  grid  point. 

(iii)  For  i  ■  1  ...  m  consider  the  intersection  of  this  plane 

with  xq  +  ^  **  0  and  use  the  resulting  linear  restriction.  For 
i  =  m  +  1  use  the  tangent  plane  itself  as  a  linear  boundary  to  the 
feasible  space. 

The  general  equation  of  a  tangent  plane  to  the  surface 
xq  +  =  f(x)  at  a  grid  point  x •  is. given  by 


n  +  1 


n 
=  i 

J-i 


3  f(x) 


a  x 


3  j 


(xJ  " 


f(x*) 


Tn.us  for  each  of  the  restrictions  f^{x)  <  0  i  «  1  ...  m  we  have 
a  set  of  linear  restrictions  of  the  form 


n  3  f . ( x) 


(Xj  -  Xj)  +  f^x")  <  0 


Similarly  for  the  restriction  x  .  .  -  g(x)  <  0  we  get  the 

n  +  -L 

set  of  linear  restrictions 


n  +  1 


n 
-  I 


a_six) 

a  x. 


] 
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x*)  -  g(x#)  <_  0 
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At  this  point  the  additional  assumption  is  made  that  the  original 
feasible  region  be  bounded.  Thus  there  is  sor>-  value  D  such  that  the 
cube  |x.|  <  D  ,  J  «  1  ...  a  contains  the  feasible  region.  Thus  if  ve  only 

w 

allow  grid  points  in  this  cube  the  linear  problem  formulated  above  hac  a 
finite  number  of  restrictions  and  the  simplex  method  could,  at  least 
theoretically  be  applied,  (in  practice  it  has  not  been  found  necessary  to 
specify  D.) 

The  concavity  of  g(x)  on  a  bounded  space  allows  us  to  postulate  a 
value  M  such  that  the  optimum  value  of  g(x)  <_  M. 

The  linear  problem  formulated  above  can  now  be  summarized  as  follows: 
max  d'x 
s.t. 

A  X  <  C 

x.  >  0  i  »  1  . . .  n 

i  - 

x  .  ,  unrestrained 
n  ♦  l 

Here  *  (x.  x,  ...  x  .  *.  ,  ,)  and  d'  *  {0  0  ...  01)  are  n  +  1 
12  n  n  ♦  l 

vectors  and  A  and  c  c an  be  written  in  partitioned  form  as  follows: 


*1  0 " 
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For  i  *  1  ...  in.  A.  is  a  matrix  whose  rows  are  formed  by  evaluating  the 

/» 1  *  t  \ 

vector  ^  ,  ...  | -  j  at  all  grid  points.  Thus  A^  is  N  x  n 

1  n  / 

where  N  is  the  number  of  grid  points.  Corresponding  to  we  have 


the  N  -  vector  Cj  whose  elements  are  given  by 


n  3 

t  y-i  »  -  f,(»*) 

4-1  ^4j..  3  1 


evaluated  at  the  grid  points  x*. 

Similarly,  A  ^  is  defined  by  the  rows 

ZQ  ♦  X 


(-H; 


» *„ / 


and  c  .  is  its  corresponding  vector  with  elements 

1 ft  +  1 


g(x») 


n  .  -i 

r  ff- 

HI  XjJ 


The  (n  ♦  l)  st  column  of  A  contains  zeros  in  all  row3  except  those 

corresponding  to  A  ,  ,  and  the  last  row  which  contains  ones.  Hote  that 
m  +  i 

this  last  row  corresponds  to  the  restriction  xn  +  ^  £  M. 

The  Dual  Problem  and  Its  Simplex  Tableau 
The  dual  to  the  linear  problem  described  above  is  given  by 
minimize  c'w 
s.t. 

A'w  d 
w  >_  0 

Here  for  convenience  the  restrictions  have  all  been  written  as  A'w  d  but 

in  fact  since  x  ,  is  unrestrained  the  last  restriction  is  an  equality, 
n  ♦  1 


,  jditKJui:  ■ 
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Thls  restriction  has  coefficients  Tillich  are  the  last  column  of  A  and  has 
a  one  on  the  right  hand  side. 

This  dual  problem  is  best,  summarized  in  tableau  form  as  shown  belcw.  This 
tableau  gives  the  rules  for  constructing  any  desired  column.  Thus  to  generate  a 
column  of  the  type  indicated  by  columns  1  through  m+1  the  following  steps 
are  neceasary: 

(a)  Rows  1  through  n  are  obtained  by  evaluating  the  indicated  partial 
derivatives  at  the  desired  grid  point. 

(b)  Row  0  is  obtained  by  evaluating  as  defined  above  at  the  desired 
grid  point  and  changing  its  sign. 

(c)  Row  n+1  is  0  for  i*>l  ...  m  and  1  for 
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The  Algorithm 

Either  the  original  or  the  dual  linear  problems  formulated  above  can 
theoretically  be  solved  by  the  simplex  method  but  it  is  clear  that  for  even 
moderately  large  problems  and  reasonable  grid  spacing  the  number  of  restrictions 
in  the  primal  (or  number  of  variables  in  fie  dual)  would  be  so  large  as  to 
exceed  the  capacity  of  most  computers. 

In  this  section  a  computational  procedure  is  developed  which  solves  the 
problem  by  a  simplex-like  algorithm  which  requires  the  formation  of  only  a  small 
number  of  the  tangent  planes  described  above. 

Since  the  algorithm  is  an  iterative  process ,  it  is  only  necessary  to  describe 
how  to  get  started,  how  to  proceed  from  one  step  to  the  next  and  when  to  stop. 

This  will  now  be  done  using  the  standard  terminology  and  notation  from  the 
simplex  method. 

At  any  iteration,  say  the  Xth,  a  basis  matrix  call  it  consisting  of 
a  *  2  columns  from  the  tableau  is  required,  or  more  precisely  ita  inverse 
is  required.  To  get  started  we  have  the  initial  basis  matrix  BQ  given  in 
partitioned  form  by 


The  inverse  of  BQ  is  easily  obtained  by  changing  the  sign  on  M. 

Assume  now  that  stage  k  of  the  iteration  has  been  reached  and  is 

available.  Denote  the  elements  of  by  b1^  ,  ij  «  0,  ...  n  +  1.  Recall  that 

the  first  row  of  B  is  the  current  pricing  vector  end  that  b^°  •  1  hence  the 

o 

current  pricing  vector  la  (1,  ...  b^D  •  Now  assuming  for  the  moment 

that  ell  columns  in  the  tableau  have  been  formed,  the  simplex  method  rev>ires  Ibe 
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computation  of  the  inner  product  of  the  pricing  vector  with  a  column  of  the 
tableau  to  see  if  this  colunn  is  eligible  to  come  into  the  basis.  Looking  first 
at  columns  of  the  type  i  »  1  ...  m  this  computation  yields  the  net  price 


* 


♦ 


£  b0j 
>1 


for  a  typical  tangent  plane  art  x*.  If  this  net  price  is  0  this  vector  is 
eligible  to  cone  into  the  basis.  Since  it  is  inconceivable  to  have  all  such 
vectors  available  for  inspection  the  following  question  is  raised: 

i.U 

Are  there  any  eligible  tangent  planes  representing  the  restriction 

and  if  so  can  the  one  with  maximum  net  price  be  identified  and  constructed? 

Looking  at  the  above  net  price  as  a  function  of  the  unknown  grid  point  and 

a 

recalling  the  definitions  of  ve  see  that  the  net  price  for  any  tangent  plane 
to  f^x)  is  given  by 


n  3  f 
NP(  x* )  *  E  r — S 

j-i  S:t.jj 


(b0'5  -  X*)  *  fj(x») 


Thus  the  net  price  is  Just  the  tangent  plane  formed  from  grid  point  x#  evaluated 
at  Xj  *  h0"^.  Now  among  all  grid  points,  this  ordinate  is  maximum  for  x*  =  b0*^ 
and  the  magnitude  of  the  net  price  is  f\(b0^).  Thus  among  all  the  vectors  which 
can  he  formed  from  f^Cx),  the  one  with  maximum  net  price  is  formed  by  allowing 
the  grid  point  to  have  coordinates  x^  =  and  further,  this  maximum  net  price 
is  given  by  evaluating  f^  at  this  grid  point.  Thus  if  f^b0^)  is  negative  there 
are  currently  no  columns  in  the  ith  set  which  can  come  in  and  another  set  must 
be  inspected. 

Similarly,  if  the  pricing  vector  is  applied  to  the  (m  ♦  l)  et  set  of 
restrictions  the  net  price  is 

>  n 


a  ♦  1 


I  ffi- 
J-l  X3  . 


3°J  ♦  b°  n  +  1 
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which  upon  rearrangement  can  be  written  as 

»(*.)  .  b"  ^(b«J  -  xj,  ♦ 

VtllCh  1S  •*ust  b°  M1  -  <th'  oratoKe  «  Xj  -  b0J  to  tht  tregent  plane  fon»d 
from  grid  point  **) .  Again  this  net  price  is  maximum  for  x*  -  b0J  and  the 
maximum  net  price  is  given  by  b°  n+1  -  g(b). 

aus  in  summary,  the  pricing  operation  consists  simply  of  evaluating  the 
functions  f^x),  U1  ...  m  and  xn+1  -  g(x)  at  Xj  .  b°J,  j-l  . . .  n+1.  If  any 
of  these  is  positive  then  the  corresponding  vector  can  be  brought  in.  In 
practice  the  one  corresponding  to  the  largest  net  price  might  be  selected.  The 
vector  is  then  formed,  up-dated  and  brought  into  the  basis  according  to  the 
usual  simplex  procedure  to  get  B^  and  the  process  is  repeated.  The  steps 
are  repeated  until  all  net  prices  are  <  0,  or  more  reasonably,  until  no  net 
price  is  >  e  where  e  is  some  small  positive  number. 

According  to  the  dual  theorem  of  linear  programming  the  pricing  vector  for 
the  terminal  B  1  matrix  gives  the  optimum  solution.  That  is,  x.  =  b°J  and 
g(x)  =  b°  n+1.  J 


Observations 

The  sequence  of  pricing  vectors  generated  as  described  above  is,  in  fact,  a 
sequence  of  trial  solutions  which  converge  on  the  true  solution  from  outside  the 
feasible  region.  Each  trial  solution  is  at  the  intersection  of  the  n+1 
tangent  planes  which  are  currently  in  the  basis. 

The  choice  of  value  of  M  requires  seme  knowledge  of  the  problem  although 
the  previous  paragraph  shows  that  if  M  is  too  small  this  will  be  indicated  by 

the  fact  that  b°  n+1  =  M  at  termination.  There  is  theoretically  nothing 
wrong  with  choosing  M  too  large. 
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Some  problems  have  been  encountered  in  which  some  of  the  hypotheses  nuve 
not  been  satisfied  in  certain  regions.  For  exeng>le,  /" x™  ia  concave  for 

0  but  not  real-valued  for  x^  <  0.  Thus  if  the  pricing  vector  at  some 
stage  has  b°*  <  0  then  the  pricing  operation  can  not  be  performed.  However,  in 
this  case,  the  restriction  ve<rfcor  in  the  tableftu  vould  be 

eligible  to  cone  in  and  if  brought  in,  will  remove  the  difficulty. 

In  general,  if  ail  assumptions  are  not  satisfied  it  may  be  possible  to 
apply  the  algorithm  tut  care  should  be  taken. 
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PART  II:  Documentation  of  „he  Pro,-ra.o 

I.  Description  of  the  Problem 

As  indicated  In. JPert  I  the  problem  to  be  solved  is: 

maximize  g(x) 

subject  to  ft(x)  <0,  i  =  1,  m  , 

where  g(x)  is  concave  and  the  f^(x)  are  convex,  real  valued  functions  of  the 
n-vector  x  for  all  real  x,  and  the  functions  are  differentiable.  Mote  that 
description  as  a  minimization  problem  is  also  efforded,  since  minimizing  -g(x)  is 
equivalent  to  maximizing  g(x). 

For  simplicity  of  programming  we  assume  that  the  constraints  are  segmented 
into  four  categories  in  the  following  order: 

lower  bound  constraints  denoted  by  Xj^  >  1^,  i  =  1,  NLB 

upper  bound  constraints  denoted  by  x^  <  u. ,  1  =  1,  ...»  HUB 

other  linear  constraints  denoted  by  Dx  <b,  D  is  (NLIN)  x  (NVAR) 

other  convex  constraints  denoted  by  hi  <  0,  i  =  1 . NCOK 

Letting  the  functional  be  denoted  by  h^^  +  ^  the  problem  can  be  described 
as 


“ICl"1'e  VOH  .  1<«> 

subject  to 

-xi  +  \  1  o. 

i«  ■■  •  «  *  y  NLB 

(i) 

xi  -  uL  <  0, 

i  *  1)  * • * )  NUB 

(2) 

Dx  -  b  <0, 

where  D  ts  (N  LIN  x  N  VAR) 

(3) 

and  b  is  (NLIN  x  l) 
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h1(x)  <0,  1*1,...,  NCOS 


H.  Description  of  the  Program 


CONVEX  Is  a  subprogram  written  in  FORTRAN  IV  (G)  which  executes  the 
computational  procedure  consisting  cf  simplex-like  iterations 
with  special  pricing  operations.  The  user  is  required  to  write  a  main  program 
and  function  subprogram  which,  together,  specify  the  particular  problem  CONVEX 
is  to  solve,  and  call  CONVEX. 

The  main  program  will  have  available  to  CONVEX  fa  array  consisting  of  the 
lower  and  upper  bounds  and  the  left  side  of  (3).  This  information  on  upper-and 
lower  bounds  will  be  provided  through  the  single  dimensioned  arrays  ILB,  ELB, 
IUB,  BUB,  where 


ILB  gives  tho  subscripts  of  the  variables  having  lower  bounds 
BLB  gives  the  lower  bound  for  the  ILBtb  variable 
IUB  gives  the  subscripts  of  the  variables  having  upper  bounds 
BUB  gives  the  upper  bound  foi  the  FJBth  variable. 


For  example,  consider  a  seven-variable  problem  having  upper  and  lower  bounds 
as  follows : 


^X3^3 

<  U^  or,  if.  proper  input  fora. 


"*3  +  1*3  —  O 


-x6  +  l6  <  0 


^6  -  *6 


X3-U3<0 


x5  -  U5  <  0 


This  information  may  be  provided  by  having 


I 
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hb(i)  >  3 

ILB(2)  =  6 
BLB(l)  -  Lj 
BLB(2)  =  Lg 
IUB(1)  =  3 
IUB(c)  =  5 
BUB(l)  =  -U3 
BUB  (2)  =  -U? 

The  other  linear  constraints  are  entered  as 


a(i,  j)  =  [d  :  -bT, 

* 

(  x  \ 

vfaere  Ai-j-J  <  0  represents  constraints  as  in  (3).  Values  of  the  following 
variables  should  also  be  supplied: 


NVAR:  number  of  variables  . 

NET:  total  number  of  constraints  plus  one  • 

MLB:  number  of  lower  bound  constraints. 

NUB:  number  of  upper  bound  constraints. 

NLIN:  number  of  other  linear  constraints  • 

NCCNs  number  of  non-linear  constraints  • 

BICW:  upper  bound  on  objective  function. 

DELTAM:  amount  by  which  BIGM  is  increased  if  the  original  estimate  of 

31(34  is  too  small. 

NALTER:  If  NALTEh  =  0,  alternate  net  pricing  is  employed.  Otherwise, 

albernate  net  pricing  is  employed  except  that  priority  restriction, 
are  checked  each  iteration.  If  NALTER  -  1,  "priority"  tteer.s 
lower  bounds.  If  NALTER  =  2,  "priority"  means  lower  and  utper 


EPS1: 


EPS2: 


KEPS2: 


EPSPIV: 


IBAS: 


BIHV  t 


X: 


DELTAX: 


I MBAS: 
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bounds.  If  IIA1TER  "  3,  "priority"  means  lower  and  upper  bounds 
and  other  linear  constraints.  If  halter  .  4,  "priority"  means 
all  constraints. 

represents  the  stopping  criterion  associated  vith  the  net  pricing 

operation,  i.e.,  the  program  may  stop  if  the  maximum  net  price 
is  less  than  EPS  1. 

represents  the  stopping  criterion  associated  with  the  sum  of  the 

absolute  value  of  changes  in  the  values  of  the  unknowns  from 

NVAfi 

iteration  t to  iteration t+k,  i.e.,  z  |xb  -  x*+k|  <E?S2. 
represents  k  in  the  description  of  the  stopping  criterion 
associated  with  EPS2,  and  must  be  less  than  or  equal  to  15. 
represents  pivot  check  stopping  criterion,  i.e.,  the  program 

may  stop  if  maximum  element  in  entering  vector  is  less  than 
EPSPIV. 

Instructs  program  concerning  which  basis  to  expect, 

and  must  be  0,  1,  2,  or  3-  For  IBAS  =  0,  the  initial  basis  inverse 

is  BIHV,  supplied  by  the  user.  For  IBAS  =  1,  the  lower  bound 

constraints  are  automatically  used.  For  IBAS  =2,  the  upper 

bound  constraints  are  automatically  used.  For  IBAS  =  3,  the 

co  ordinates  of  an  initial  feasible  point  are  supplied  in  X. 

used  only  if  IBAS  =  0,  in  which  case  Birr/  contains  the  inverse 

of  the  initial  basis  and  has  dimension  (NVAR  +  1)  x  (NVAR  +  l). 

used  only  if  IBAS  =  0  or  IBAS  =  3,  in  which  case  X  contains  the 

current  solution  vector  but  has  dimension  NVAR. 

used  only  if  IBAS  *  3,  in  which  case  it  is  the  amount  by  which 

any  artificial  bound  is  incremented  if  it  appears  in  a  terminal 
solution. 

used  only 

the  index 


if  IBAS  =  0,  in  which  case  the  user  must  supply  in  UlBAS 
of  the  constraint  represented  in  the  Tth  .... 
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INVERT : 


ITFRIiT : 


JSIH: 


instructs  program  as  to  hew  many  Iterations  to  perform  before 
re -inverting. 

^eHs  program  how  many  iterations  are  to  be  performed  before 
intermediate  results  are  printed. 

If  JBIN  =»  0,  B  1  is  included  in  the  information  printed  out. 
Otherwise ,  it  is  excluded. 


In  addition,  the  main  program  must  include  the  following  three  statements: 

DIMENSION  BINV(60,  60),  X(60),  INBAS(6o) 

COi\J£>N/BLK  i/a(6o,  6o),  blb(6o),  bjb(6o),  ilb(6o),  iub(6o) 
comioh/blk  2/nvar,  nlb,  nub,  nun,  ncon,  net  * 

A  function  subprogram  of  the  form  SETUF(l,  J,  X,  C)  is  required  to  provide 

V  1  "  lf  •*’»  *!0H  +  1*  and  a|“hi»  1  *  •••>  N20N  +  1;  J  =  1,  NVAS.  For 

specified  I,  SETUP  will  provide  h^x)}^  when  J  -  0,  and  g|-  h^x)^ 

J  s  1>  •••>  NVAR  when  J  a  1,  ^ 

III .  An  Example 

Consider  the  following  problem: 

Maximize  +  Xg 

subject  to  -,7X1  +  Xg  <  1 


2X1  +  ^4  -  6 

0  <  Xg  <  1.3 


x1  >0 


*For  the  double-precision  veraion,  these  arrays  should  be  dimensioned  30,  or 
30  x  30,  except  for  As  which  should  be  50  x  30.  * 
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The  main  program  may  appear  as  follows; 


& 


DIMENSION  31157(60,  60),  x(6o),  INBAS(60) 

cotaJDN/sLK  i/a(6o,6o)  ,  813(60),  bub(6o),  &5(6o),  IUn(50) 

C0M40N/BLK  2 /nvar,  nlb,  nub,  nlin,  noon,  net 

WAR  =  2 

NET  =  6 

NLB  =  2 

NUB  *  1 
NLIN  =  1 
NCON  =  1 
BIGM  =  4. 

DELTAM  =  1. 

NALTER  -  0 
EPS1  =  .001 
EPS2  *  .001 
EPSPIV  =  .0001 
INVERT  =  50 


NEPS2  =  10 


-l6f 


TEAS  =  1 
ITPRRT  =>  1 

vTEIN  =  0 

ILB(l)  ■  1 

ILB(2)  -  2 

BIB(l)  »  0. 

BLB(2)  *  0. 

IUB(l)  »  2 

B0B(1)  =  -1.3 

A(l,  1)  =  -.7 
A(l,  2)  »  1. 

A(l,  3)  »  -1. 

CALL  CONVEX  (IBAS,  BINV,  X,  TNBAS,  NALTER,  EPS1,  EPS 2,  NEPS2,  EPSPIV, 
^INVERT,  ITFRNT,  JBIH,  BIGM,  DELTAM,  DELTAX) 

STOP 

END 

SETUP  may  appear  as  follows: 

FUNCTION  SETUP  (I,  J,  X,  COL) 

DIMENSION  X(6o) ,  C0L(60) 

SETUP  =  0. 

IF  (J.  EQ.  0)  GO  TO  (100,  200),  I 

IF  (J.  EQ.  1)  GO  TO  (300,  400),  I 

100  SETUP  *  2.*X(l)**2  +  3-*X(2)**2  -  6. 

RETURN 

200  SETUP  -  2.*X(l)  -  X(l)**2  +  X(2) 


RETURN 


300  C0L(1)  -  4.  *  X(l) 

C0L(2)  *  6.  #X{2) 

RETURN 

400  COL(l)  =2.-2.  »£l) 

C0Xi(2)  =  1  • 

RETURN 

END 

The  program  cat.  then  be  assembled  as, 
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Institute  of  Statistics 
Texas  A&M  University 
College  Station,  Texas  77^3 


% 


Sg- 

I 


s 

I 


Mr.  Ernest  Davis 
Department  of  Ag.  Economics 
Texas  A&M  University 
College  Station,  Texas  77$+3 


s 

f-' 

a 


KSY  XOHO.I 


LIMK  • 


1 


Convex 


Progr  aiming 
Optimization 


Onclasalfled 


Security 


)D  .'"“..1473 


«»*  (i  I 

PAGE  2) 


( BACK ) 


